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Abstract 

We consider two simple models of higher-derivative and nonlocal quantum 
systems. It is shown that, contrary to some claims found in literature, they can 
be made unitary 
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Considerable attention has been paid recently to the problem of quantizing the 
theories which contain higher derivatives or are even nonlocal in time. This is mainly 
due to the fact that such theories arise naturally when one attempts to model the 
dynamics on noncommutative space-time with the help of star product construction [1] 
-T- [3]. It has been shown that the field theories with space-time noncommutativity lead 
to nonunitarity and acausality, at least when quantized with the help of naive Feynman 
rules [I] -T- [6]. There are several proposals of alternative quantization schemes which 
seem to cure the trouble with nonunitarity [7j -j- [10] . However, they give rise to 
other problems [TT] , p2] so the question of the existence of consistent quantum theory 
remains open. In view of that one is tempted to study simple models to gain more 
understanding of the problems we are faced with. However, even in the simplest cases 
there is some misunderstanding concerning the consistency of the relevant models. In 
order to clarify some issues we study two simple systems which are sometimes claimed 
to lead to nonunitary evolution. 

Let us first consider the Lagrangian 



L = \(i?-^) (1) 

The corresponding classical equation of motion reads 

q (iV) _tf q = ( 2 ) 

The general solution to eq.pi) is a linear combination of exp(±iflt) and exp(±Qt). It 
is sometimes claimed [13] that due to the existence of exponentially rising solutions 
exp(±Qt) the quantum Hamiltonian of such type of systems is not hermitean and the 
evolution operator exp (— f^-) is not unitary. 

Let us consider this problem in some detail. In order to quantize the model we have 
to put it first into Hamiltonian form. This can be done with the use of Ostrogradsky 
formalism [T3] [16]. One defines the canonical variables 

qx = q, q 2 = q 

5L 8L d (dL\ 

U ^Tq=^~di{^)=- q (3) 

= 6L _ dL _ .. 

2 8q dq 

and the Hamiltonian 

H = U iqi + U 2 q 2 — L — IL iq2 + ^Tlj + ^q\ (4) 

It is straightforward to check that the canonical equations of motion imply the initial 
equation for q. 

Our system, being linear, can be immediately quantized. The question we intend to 
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solve is whether the resulting quantum system is unitary. The argument against uni- 
tarity based on the existence of exponentially growing solutions to classical equations 
of motion might go as follows. Let 

X = Ui - im 2 - Sfqi - tffc (5) 

then 

[x, h] = ihnx (6) 

(which is equivalent to the existence of exponentially growing solutions to Heisenberg 
equations of motion). Assume now that H is selfadjoint and let | E) be the eigenvectors 
of H (possibly the generalized ones, corresponding to continous spectrum). Then eq. 
([6]) implies 

(E- E' - ihtt) (E' | X | E) = (7) 
Eq. ([7]) is fulfiled only provided 

(E' | X | E) = (8) 

for all E, E' which in turn means that the canonical variables are lineary dependent 
contradicting the canonical commutation rules. 

The only weak point in the above reasoning is the assumption on the existence 
(even in distributional sense) of the matrix elements (E' \ X | E). In fact, we shall 
show below that the theory is unitary while (E' \ X \ E) are not well defined. 
To this end we make a simple canonical (hence unitary) transformation 



Ul = ^\ Pl + P2 )' n 2 = -j=(-x 1 + x 2 ) (9) 

Then the Hamiltonian takes the following form 

Ms* + -G*- £*0 (10) 

H is the difference of harmonic oscillator and the inverted harmonic oscillator. Both 
terms depend on different variables so if they are selfadjoint and generate unitary 
dynamics so is H. Only the second piece calls for some explanation. It is unbounded 
from below but only quadratically. Therefore, it defines selfadjoint operator with purely 
continuous spectrum covering the whole real axis |17j . 

In order to find the coordinate representation for eigenvectors of H we have only 
to find the corresponding wave functions for inverted harmonic oscilator. This is quite 
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straightforward [18] and results in the following final formula for eigenfunctions of H. 

(Xt,X 2 | n,e;±) = (11) 



y/ir2 n n\ 

V h ^ch{^) m ~ 2 \ V n \ v/2 J 

here H n are Hermite polynomials while D u - the parabolic cylinder functions. The 
energies of the above states are given by 

E = E(n,e) = hn(n + - S \-e (12) 

Consider now the matrix elements (E' | X | E). By virtue of eqs. (jSJ) and (Q we find 

X = V2VL(P 2 - ttX 2 ) (13) 

Now, using ( fTTj) . ( fTBl and the asymptotic form of parabolic cylinder functions [18] we 
find that the matrix elements (n', e' , X' | X | n, e, A), A, A' = ±1, are expressed, apart 
from regular contributions, in terms of badly divergent integrals of the form 

dxe^ 1 ™ ~ j dye'i 1 ^) (14) 

Therefore, the matrix elements under consideration are not well defined even in the 
distributional sense. 

As a second model we consider the nonlocal harmonic oscillator 
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Euler-Lagrange equation 



L = U(t) 2 -^-q(t-^)q(t+ ^) (15) 



& S -TTT = ( 16 ) 

Sq(t) 1 ; 

reads 

q{t) + ^-{q{t-T)+q{t + T)) = U (17) 

This theory can be quantized using the method of Pais and Uhlenbeck [21J. To this 
end we rewrite the Lagrangian in equivalent form 

L = -\q(t)q(t) - ^((q(t)q(t — T) + q(t)q(t + T)) (18) 



or 



1 , x / ( d\ 2 9 , / „ d 



£ =-2*'(UJ Wc T*Jr (<) (19) 

According to the prescription given in Ref. [2T] we consider the entire function 

$0) = z 2 + uj 2 ch(Tz), u = z 2 (20) 

It is not difficult to verify that: (i) $(«) has no multiple zeros except a discrete set 
of values of uT; in the latter case the double zeros are real and negative, (ii) there is 
a finite nonempty set of real negative zeros, z 2 = —£l 2 , i = 1, ...,m, m > 1, (m) 
there is an infinite number of complex pairwise conjugated zeros, z 2 = —uj 2 , z 2 = 

— ZD 2 ., k = 1,2, ... and u k * s absolutely convergent. 

fc=i 

By virtue of the above properties one can write 
Following [2T] we define new variables 

ft 3 n (i + ^) n (i + ^) (i + « 

™ / 1 \ ^ / 1 d 2 ^/ 1 d 2 \ 

* S S( 1+ ^ S( 1+ ^* S( 1+ e?* 9 (22) 

With the above definitions the Lagrangian takes the form 

1 m ~ / d 2 \ ~ 

L =-2£^<Ua+n?JQi- (23) 

By rescaling Qi — >■ 9*— , Qfc — >■ Qa. — > -%=, and passing to the Hamiltonian 

formalism we find 

-i m i oo 

h = 5 ^ 2 + ^ 2 ) + o E((^ 2 + ^ 2 Q 2 ) + (Pi + (24) 

z i=i z fe=i 
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This is, however, not the end of the story because the variables Q k ,P k are complex. 
In order to find the relevant real variables we perform the following complex canonical 
transformation [21] 

Pk = i^V^k ((pife + iqik) + i(P2k - z<?2fc)) 

Pk = ^V^k{{Pik ~ iqik) ~ i{P2k + iq2k)) (25) 

1 

Qk = 7T^= ((Pi* ~ ^lfc) + i(P2k + iq2k)) 



— 2 

Qk = 7T7^ ((pi* + iqik) - i{P2k - iq2k)) 

which allows to rewrite eq. ff24l) as 

H = (26) 

\ m _ 00 

= o + ^Q 2 i) - J2(Pniu k (q lk pi k + q 2k p2k) + Reuj k (q lk p 2k - q2kP\k)) 

^ i=l k=l 

So is a sum (with alternating sign) of finite number of harmonic oscillators and 
an infinite number of terms which are linear combinations of dilatation and angular 
momentum in two dimensions. All terms depend on different variables and can be 
easily quantized; only dilatation calls for ordering rule - we adopt the simplest one: 
qiPi + <M?2 -> + Ml + hP2 + P2&)- 

Let us consider one term of the second sum on the RHS of eq. fl2"61) . It is of the form 

h = ^(<?iPi + Pi 9i + + P2Q2) + K<7iP2 - 92P1) (27) 

Dilatations commute with angular momentum so the spectrum of h reads 

£n,\ = fJ>X+ vn, n — 0, ±1, —00 < A < 00 (28) 

and the corresponding wave functions (for example, in coordinate representation) can 
be easily found [2"T| . 

Now, let us consider classical equations of motion implied by h. They read 

ft = /Uft - ve ik q k , pi = -fxpi + ve ik p k (29) 

The solutions are linear combinations of exp((/j, ± iu)t) (for </s ) or exp((—/j, ± 
(for p's ). Therefore, we are dealing with exponentially growing solutions. Following 
the same line of reasoning as in the first example one could conclude that the energy 
must take complex values. Again this conclusion results from the assumption that the 
matrix elements of canonical variables in the energy representation are well-defined 
while it is straightforward to check that they are badly divergent (cf. eq.(34) of Ref. 



Let us now make few remarks about path integral approach. Some care must be 
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excercised when dealing with the Hamiltonians unbounded from below. Let us take 
as an example the inverted harmonic oscilator with the "frequency" u. As we have 
mentioned above the Hamiltonian is here a well-defined selfadjoint operator with purely 
continuous spectrum extending over the whole real line. 
The propagator function 

oo 

K(x,y;t)=Y, J d^—y^Ex(xW E x(y) (30) 

^=±1—00 

cannot be continued to imaginary time whatever the sign of t is. On the level of path 
integration the Euclidean path integral is not well-defined due to the fact that the 
integrand ceases to be of the form of the exponent of negative definite functional. 
However, we can do the path integral directly (without refering to imaginary time 
approach) by applying Trotter formula and doing Fresnel integrals. The result, when 
extended to whole time axis, reads 



K(z,y;t) = (1 ~ iS9nt) J ^e^tr^+y 2 ^- 2 ^ (31) 
2jnh 1*1 V shut 



This result can be checked as follows. By virtue of eq. ( 1301) one gets 

2tt 



/ dte^K(0,0 ]t ) = - £ |^ A (0) | 2 (32) 

4c U A=±l 

Using eq. (I3T]) and the explicit form of *f?§x ( c ^- ec L- CCQ) ) we find that (1321) holds true. 
Consider now K(0, 0; t) for t > 0: 

Let us forget for a moment that the imaginary time representation is ill-defined. In 
order to extract the information about the spectrum of our Hamiltonian we continue 
( |33i) to imaginary time and consider the t — > oo behaviour. We find that K(0,0;t) 
becomes periodic in Euclidean regime which means that the energies are imaginary! 

For our first example the propagator is simply the product of two propagators, one 
for harmonic oscillator and one for inverted harmonic oscillator. Therefore, we see 
that the inproper application of path integral method will result in wrong conclusions 
concerning the energy spectrum and unitarity. 

The same applies to the nonlocal harmonic oscillator. We found that the Hamil- 
tonian is a sum of commuting pieces, many of them having the spectrum unbounded 
from below. Again, it appears that an attempt to derive the energy eigenvalues from 
Euclidean version of path integral leads to complex energies, i.e. to the conclusion that 
the theory is not unitary. 
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